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The lifetime of two-vibron bound states in the overtone re- 
gion of a one- dimensional anharmonic molecular lattice is in- 
vestigated. The anharmonicity, introduced within an attrac- 
tive Hubbard Hamiltonian for bosons, is responsible for the 
formation of bound states which belong to a finite linewidth 
band located below the continuum of two-vibron free states. 
The decay of these bound states into either bound or free 
states, is described by considering the coupling between the 
vibrons and a thermal bath formed by a set of low frequency 
acoustic phonons. The relaxation rate is expressed in terms of 
the spectral distribution of the vibron/phonon coupling and 
of the two-vibron Green operator which is calculated exactly 
by using the number states method. The behavior of the two- 
vibron bound states relaxation rate is analyzed with a special 
emphasis on the influence of the anharmonicity. It is shown 
that the rate exhibits two distinct regimes depending on the 
thermal bath dimension. When the bath dimension is equal 
to unity, the rate increases with the anharmonicity and the 
decay of the two-vibron bound states into the other bound 
states appears as the main contribution of the rate. By con- 
trast, when the bath dimension is equal to 2 and 3, the rate 
decreases as the anharmonicity increases indicating that the 
two-vibron bound states decay into the two-vibron free states 
continuum. 



I. INTRODUCTION 

The nonlinear nature of vibrational excitations in low- 
dimensional molecular lattices plays a key role for energy 
transfer as well as energy storage in both physical, chem- 
ical and biological systems. 

This feature was first pointed out by Davydov and co- 
workers [1] to explain how the energy released by the 
hydrolysis of ATP, partially stored in the high frequency 
C=0 vibration (amide-I) of a peptide group of a protein, 
can be transported from one end of the protein to the 
other. The dipole-dipole coupling between the different 
peptide groups leads to the derealization of the internal 
vibrations and to the formation of vibrons. The interac- 
tion between the vibrons and the phonons of the protein 
induces a nonlinear dynamics which counterbalances the 
dispersion created by the dipole-dipole coupling. It yields 
the creation of the so called Davydov's soliton which ap- 
pears as the solution of the Nonlinear Schrodinger (NLS) 
equation within the continuum approximation. Never- 
theless, it has been shown that several properties, such 
as the space dimension [2] and the discreteness of the 



lattice [3], tend to destroy the integrability of the NLS 
equation. 

However, a remarkable feature given by the discrete 
version of this equation is the occurrence of intrinsic 
localized modes or discrete breathers (for a recent re- 
view, see for instance Ref. [4]). Contrary to solitons, 
intrinsic localized modes do not require integrability for 
their existence and stability. They are not restricted to 
one-dimensional lattices and it has been suggested that 
they should correspond to quite general and robust so- 
lutions [5]. In molecular lattices, discrete breathers can 
be induced by the intramolecular anharmonicity of the 
molecules. Indeed, within the semi-classical approxima- 
tion [6] , the vibron dynamics is described by the discrete 
NLS equation. As a result, the presence of an intrinsic 
localized mode corresponds to a local accumulation of 
the internal energy, which might be pinned in the lattice 
or may travel through it. Breather excitations are thus 
expected to be of fundamental importance for both vibra- 
tional energy storage and transport in molecular lattices. 

In spite of the great interest that breathers have at- 
tracted, no clear evidence has yet been found for their 
existence in real molecular lattices. By contrast, two- 
vibron bound states (TVBS), which can be viewed as 
quantum breathers, have been observed in several sys- 
tems. Indeed, in quantum lattices, the intramolecular 
anharmonicity, which is responsible for the interaction 
between the vibrons, favors the formation of bound states 
[7] . When two vibrons are excited, a bound state corre- 
sponds to the trapping of two quanta over only a few 
neighboring molecules with a resulting energy which is 
lesser than the energy of two quanta lying far apart. The 
lateral interaction yields a motion of such a state from 
one molecule to another, thus leading to the formation 
of a dclocalized wave packet with a well-defined momen- 
tum. As mentioned by Scott [8], these bound states can 
thus be viewed as the quantum counterpart of breathers 
or soliton excitations. From an experimental point of 
view, the formation of TVBS in low-dimensional molec- 
ular lattices was observed in molecular adsorbates such 
asH/Si(lll) [9], H/C(lll) [10], CO/NaCl(100) [11] and 
CO/Ru(100) [12] using optical probes. Bound states in 
the system H/Ni(lll) were investigated using high reso- 
lution electron energy loss spectroscopy [13]. In addition, 
the powerfulness of two-vibron infrared spectroscopy in 
probing the structure of admonolayers has been recently 
demonstrated both experimentally [14] and theoretically 
[15]. 
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The stability of both classical and quantum breathers 
under the influence of a fluctuating thermal bath is of 
fundamental importance for energy storage and energy 
transport at finite temperature. For instance, at biologi- 
cal temperature, the lifetime of the Davydov's soliton is 
still an open question. Indeed, it has been shown that 
the amide-I excitation, in vivo, corresponds to a localized 
state [16,17]. Instead of traveling in a coherent manner, 
it follows a stochastic, diffusional path along the lattice. 
In other words, the single-vibron Davydov soliton does 
not last long enough to be useful at biological temper- 
atures and multi-vibron solitons have been proposed as 
good candidates [18,19]. In the same way, the influence of 
white noise and nonlinear damping on the stability of the 
breathers of the discrete NLS equation has been studied 
[20] . It has been shown that the lifetime of the breathers 
increases as the nonlinearity increases and is inversely 
proportional to the noise variance. By contrast, for quan- 
tum lattices, it was pointed out by Jakob et al. [21] that 
the dephasing constant of TVBS is larger than the single- 
vibron dephasing constant by a factor four. This feature, 
which occurs when both kinds of states interact in the 
same way with the thermal bath, leads to a decrease 
of the TVBS lifetime when compared with the single- 
vibron lifetime. Similar results were previously obtained 
by Bogani et al. [22] who have characterized the relax- 
ation of bound states in the combination regions of two 
internal modes in a molecular crystal. The authors have 
shown that the lifetime of a bound state decreases as the 
anharmonicity increases, excepted at very low tempera- 
ture. They assumed that the anharmonicity yields the 
occurrence of a bound state with zero linewidth and that 
the relaxation, mediated by the phonons of the crystal, 
leads to the decay of the bound state into the two-vibron 
continuum. 

In this paper, a theoretical study is performed to in- 
vestigate the TVBS lifetime in the overtone region of 
a one-dimensional molecular lattice. The bound states, 
which are described by using a local approach, belong to 
a finite linewidth energy band located below the contin- 
uum of two-vibron free states. The relaxation processes 
are considered by introducing a coupling between the vi- 
brons and a thermal bath formed by a set of low fre- 
quency phonons. Such a coupling, which is assumed to 
be responsible for dephasing mechanism only, does not al- 
low for energy relaxation. In other words, the interaction 
with the thermal bath does not modify the vibron num- 
ber but induces transitions between the different two- 
vibron eigenstates. The formalism can be applied to de- 
scribe different physical situations such as the dynamics 
of low-dimensional adsorbed nanostructures (molecular 
nanowire adsorbed on stepped surfaces) and the collec- 
tive vibration of the peptide groups in a protein (vibrons 
in a- helix polypeptide proteins). In the first situation, 
the thermal bath is formed by the phonons and the li- 
brons of the adsorbed nanowire and by the phonons of 
the substrate [23-25]. In the second situation, it is char- 
acterized by the low frequency vibrational modes of the 



protein as well as by the biological surrounding. 

The paper is organized as follows. In Sec. II, the 
Hamiltonian of the vibrons and of the thermal bath, as 
well as the coupling Hamiltonian, are described. In Sec. 
Ill, we first summcrize the number states method used 
to characterize the TVBS and the equivalence between 
the two-vibron dynamics and the motion of a single ficti- 
tious particle in a one-dimensional lattice is established. 
Then, the TVBS relaxation rate due to the coupling with 
acoustic phonons is expressed in terms of the two-vibron 
Green operator which is finally determined. In Sec. IV, a 
detailed analysis of the behavior of the rate is performed 
depending on the values taken by the relevant parameters 
of the problem (the anharmonicity, the temperature, the 
thermal bath dimension). These results are finally dis- 
cussed in Sec. V. 



II. MODEL AND HAMILTONIAN 

Let us consider a set of N molecules periodically dis- 
tributed on the sites of a one-dimensional lattice. To 
characterize the vibron dynamics of this molecular lat- 
tice, we assume that each molecule n behaves as an in- 
ternal high frequency oscillator described by the stan- 
dard creation and annihilation vibron operators and 
b n . The Hamiltonian which describes the collective dy- 
namics of the internal vibrations is expressed as (using 
the convention K = 1) 

ff„ = ^/i„ + ^$(n,n>+&„' (1) 

n n,n' 

where h n defines the internal vibration Hamiltonian of 
the nth molecule and $(n, n') represents the lateral dy- 
namical coupling between the internal coordinates at the 
quadratic order. To minimize the number of parame- 
ters, we consider interactions between nearest neighbor 
molecules and define the hopping constant <&. Note that 
long range lateral interactions may affect the two-vibron 
dynamics by inducing changes in the vibron dispersion 
curve. However, these effects are expected to be rather 
weak for a one-dimensional lattice and are disregarded 
in the present work. The internal vibration Hamiltonian 
h ni including the intramolecular anharmonicity of each 
molecule, is described according to the model of Kim- 
ball et al. [7] . In this model, the intramolecular potential 
is expanded up to the fourth order with respect to the 
internal coordinate and a unitary transformation is per- 
formed to keep the vibron-conserving terms, only. The 
resulting Hamiltonian H v is essentially a Bose version of 
the Hubbard model with attractive interactions written 
as, 

H v = ^u: G b+b n - Ab+b+b n b n + <S>(b+b n+1 + b+ +1 b n ) 

n 

(2) 
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where luq stands for the internal frequency of the nth 
molecule and A denotes the positive anharmonic param- 
eter. 

These molecules interact with their surrounding which 
acts as a thermal bath and which is responsible for the 
relaxation mechanism. To mimic the influence of the sur- 
rounding of the molecules, we assume that the thermal 
bath is formed by a set of low frequency harmonic oscil- 
lators which the Hamiltonian is expressed as 

H b = %a+ a i ( 3 ) 

where a+ and a q stand for the boson operators of the q th 
mode with frequency Q g . By neglecting the population 
relaxation, the thermal bath leads to stochastic fluctu- 
ations of the dynamical parameters (internal frequency, 
anharmonic parameter, force constant) which character- 
ize the vibron dynamics. In this work, we consider that 
the main effect of the bath is to induce a random mod- 
ulation of the internal frequency of each molecule. The 
coupling Hamiltonian between the vibrons and the ther- 
mal bath is thus written as 

A^ 6 = ^Ao;„6+^ (4) 

n 

where Aui n is an operator with respect to the bath vari- 
ables. In a general way, it depends on the bath coordi- 
nates in a very complicated manner. However, to sim- 
plify our discussion we consider the linear approximation 
in which Aw„ is expressed as 

Au>„ = Xq,nO>+ + X* q ,n a Q ( 5 ) 

Q 

The previous model has been extensively used to study 
relaxation mechanisms in condensed matter (see for in- 
stance Refs. [26-29]) and it has been shown that the ef- 
fect of the thermal bath is completely determined by the 
coupling distribution J {to) defined as 

= ^ E i i 2 - ( 6 ) 

n.q 

The theory presented below to calculate the TVBS re- 
laxation rate depends on the coupling distribution and 
can thus be applied to different kinds of thermal baths. 
In the present work, we shall restrict our analysis to an 
interaction with a set of acoustic phonons. In that case, 
the coupling distribution is determined by assuming a 
Debye model for the phonon density of states coupled 
with a deformation potential approximation [26] . There- 
fore, J {to) scales as ui a 7 where a stands for the dimension 
of the phonon lattices, and is defined as, 

_ f j Q (uj/n c r ifo><tt c 



where tt c denotes the cutoff frequency of the phonon 
bath. Note that, in the present paper, the bath dimen- 
sion is considered as a parameter which can take the val- 
ues 1, 2 and 3 in order to represent different physical sit- 
uations. For instance, for an adsorbed one-dimensional 
molecular lattice, the vibrons can be coupled preferen- 
tially with the phonons of the molecular lattice. The 
resulting bath dimension is thus equal to one. In the 
same way, for a molecular adsorbate decorating a surface 
step, the vibrons can be coupled with the step phonons 
which exhibit a one-dimensional nature. Moreover, we 
can imagine other situations such as a strong coupling 
with cither surface phonons or bulk phonons for which 
the bath dimension is equal to two or three, respectively. 

III. THEORETICAL BACKGROUND 
A. Two-vibron states 

In this section, the number states method [8] is sum- 
marized in order to characterize the two-vibron states of 
the molecular lattice. Within this method, the eigen- 
states of the Hamiltonian H v (Eq.(2)) are determined by 
taking advantage of the fact that H v commutes with the 
operator ~^2 n b^b n which counts the total number of vi- 
brational quanta. As a result, the Hilbert space E of 
the vibrational states can be written as the tensor prod- 
uct E = Eo ® Ei <g> E2 <8> ... <8> E p ..., where E p denotes 
the subspace connected to the presence of p vibrational 
quanta. Single- vibron states, connected to the subspace 
Ei of dimension N, correspond to plane waves with wave 
vectors q lying inside the first Brillouin zone of the chain 
and with eigenfrequencies u{q) = uj a + 2&cos(q). 

Two-vibron states are related to the subspace E 2 , 
which the dimension N(N + l)/2 represents the number 
of ways for distributing two indistinguishable quanta onto 
N sites. A useful basis set to generate the E 2 subspace 
can be constructed by applying two creation operators 
on the vacuum state | 0) as 

I ni,n 2 ) = [1 - 5 nin2 (l - -^=)]6+6+ | 0) (8) 

where the factor in the right hand side of Eq.(8) is a 
normalization coefficient and where the restriction n 2 > 
rii is applied to avoid counting twice the vectors | ni,n 2 ) 
due to the indistinguishable nature of the quanta. As 
a result, to solve the Schrodinger equation H v | vp) = 
w I the most general two-vibron eigenstate can be 
expanded as a linear combination of the previous basis 
vectors, as 

I *> = E *("i' n2 ) I «1)«2) (9) 

"1 ,«2 

However, by taking advantage of the lattice periodicity, 
an intermediate basis set can be determined. Indeed, 
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since the two-vibron wave function ^(ni, n 2 ) is invariant 
with respect to a translation along the lattice, it depends 
on m and m = n 2 — m, and it can be written as a Bloch 
wave, as 



* (ni,n 2 ) = E *fcM^ (ril+ri2)/2 



(10) 



The total momentum k of the two-vibron states is as- 
sociated to the motion of the center of mass of the two 
quanta. The values of k are determined by imposing pe- 
riodic boundary conditions to the wave function over a 
set of TV unit cells. The resulting wave function ^fe(m) 
refers to the degree of freedom m which characterizes the 
reduced distance between the two vibronic excitations. 
It corresponds to the projection of the eigenstate | <J/) on 
the basis set | k,m) defined as, 

I k, to) = -L e jfe( " 1+m/2) I m, m + to) (11) 

Note that the previous periodic boundary conditions lead 
to limit the values taken by the index m which cannot ex- 
ceed the half size of the periodic box and thus varies from 
to (N — l)/2. Moreover, specific boundary conditions 
are introduced depending on the parity of the number 
N [8]. However, by considering the limit N — > 00, these 
features will be disregarded in this work. 

Since the total momentum k is a good quantum num- 
ber, the Hamiltonian H v appears as block diagonal in the 
intermediate basis | k,m) (Eq.(ll)) and the Schrodingcr 
equation can be solved for each k value. For a given k 
value, the representation of the restricted Hamiltonian 
H v i~ depends on the nature of the basis vectors involved 
in. Indeed, there are two different kinds of intermedi- 
ate basis vectors either describing two quanta located on 
molecules which are far apart (to > 1) or two quanta 
located onto the same molecule (to = 0). In the first 
situation, the Schrodinger equation is expressed as 

(2w - w)*fe(m) + $ fc (*k(m + 1) + * fe (m - 1)) = 

(12) 

where m > 1 and = 2$cos(/c/2). When m = 1, the 
Schrodinger equation becomes, 

(2^o - + *k*k(2) + V2$ fe * fe (0) = (13) 

whereas when m=0, it is written as 

(2w -2A-w)* fe (0) + V2l> fc *fe(l) =0 (14) 

Eqs.(12)-(14) clearly indicate how the physics of the 
two-vibron states is related to the dynamics of a ficti- 
tious single particle moving quantum mechanically in a 
one-dimensional lattice [7]. This lattice, shown in Fig. 
la, appears as a discrete, finite size chain with a length 
which extends from m = to (N — l)/2. Within this 



equivalence, the wave function ^(to) associated to the 
relative distance between the two vibrons, can be viewed 
as the wave function of the fictitious particle. Its dy- 
namics is described by a tight-binding Hamiltonian H v k 
characterized by the self-energy 2wo located on each site 
and a hopping matrix which couples nearest neighbor 
sites. The anharmonicity is responsible for the presence 
of a local defect on the site m — which induces a redshift 
—2A of the corresponding self-energy. This defect allows 
us to discriminate between two different eigenstates. The 
eigenstates of the core of the chain correspond to plane 
waves slightly perturbed by the presence of the defect 
and are connected to the free propagation of the fictitious 
particle. In terms of two-vibron states, these eigenstates 
correspond to a derealization of the separating distance 
between the two quanta and leads to the occurrence of 
two-vibron free states (TVFS) which belong to an en- 
ergy band located around 2wo (Fig. 2a). Note that if 
there is no anharmonicity (A = 0), a TVFS with a total 
momentum k = q\ + q 2 corresponds to the presence of 
two independent vibrons with wave vector qi and q 2 and 
energy w(<?i) + u{q 2 ). 

By contrast, the presence of a defect is responsible for 
the occurrence of a localized state for the fictitious par- 
ticle. This state is associated to an exponential decrease 
of the wave function from the side to = of the chain. 
In terms of two-vibron states, it corresponds to a TVBS 
which the wave function is defined as [7,8] 



tt£(m) 



AA 



(l + ^-^S m ,>)e- m ^ (15) 



2LU -LU B (k) sj2 
where Wfl(fc) stands for the TVBS frequency written as 
uj B (k) = 2u Q - 2^ A 2 + 4$ 2 cos 2 (fc/2) (16) 

and where is the length of the bond between the two 
vibrons expressed as 



! l 2uj -Lo B (k) + 2A 

4fe ~ \ 2uj - iu B {k) - 2A 



(17) 



As shown in Fig. 2a, TVBS belong to a low frequency 
band located below the TVFS continuum and character- 
ized by the dispersion curve cj B (k). 



B. Two-vibron bound state relaxation rate 

Due to the coupling Hamiltonian AH v b (Eq.(4)), 
TVBS do not represent exact eigenstates of the whole 
system "vibron+phonon bath" . More precisely, this cou- 
pling is responsible for the occurrence of transitions be- 
tween two-vibron states mediated by the emission or the 
absorption of acoustic phonons. 

By assuming the coupling Hamiltonian AH v b as a weak 
perturbation, the rate for the transition from a TVBS 
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I to another two-vibron state | (bound a = B or 
free a = F states) with frequency LO a {k') can be expressed 
by using the Fermi Golden Rule as 

W(* B -*£,) = 27r£P Q | (1rf,a | | tt£ ,0) | 2 

5(w B (fc) + fia-w*(fc')-tyj) (18) 

Eq.(18) describes a transition in the course of which the 
thermal bath evolves from an initial state | a) with fre- 
quency Cl a to a final state | (3) with frequency O/3. Since 
the thermal bath is assumed to be in thermal equilib- 
rium at temperature T, a statistics is realized over the 
initial state with the probability occupation P a and a 
sum over all the possible final bath states is performed. 
By inserting the expression of the coupling Hamiltonian 
AHyb (Eq.(4)), the total rate for leaving the bound state 
I ^ B ) obtained by summing over all possible transitions 
is expressed as 

poo 

Wgl = 2Re V / dt{Aw n (t)Aw n /(0)) 

n,n> J ° 

(*£ \b+(t)b n (t)b+(0)b n ,(0)\* B ) (19) 

where (...) stands for an average over the thermal bath 
and where the operators depend on time t according to an 
Hcisenberg representation with respect to the free Hamil- 
tonian Hq = H v + Hb. At this step, the relaxation rate 
can be written in an improved form by inserting the inter- 
mediate basis set (Eq.(ll)) for which the matrix elements 
of the population operator are defined as 

(k,m I b+b n I k',m') = ^ mm ,e^- k '> C o S {^-m) 

(20) 

Therefore, by inserting Eq.(20) into Eq.(19), the relax- 
ation rate is written as 

n.n' m./m' k' 

(Auo n {t)Auj n ,(Q))e^ k '^ n - n '\k',m | e~ iH ^ \ k',m') 
^> B *{m)^ B (m')cos((k - k')m/2)cos{{k - k')m'/2) (21) 

The TVBS relaxation rate depends on the correlation 
functions of the coupling Alu. The characteristic time of 
this rate is the correlation time r c of the thermal bath 
which corresponds to the time for which the correlation 
functions vanish. In a general way, r c is about 1 ps for 
phonons in low-dimensional molecular lattices. For in- 
stance, molecular dynamics simulations carried out for 
the system CO/NaCl found a correlation time value equal 
to 0.6 ps at 25 K [24]. We thus assume that this time 
scale is small compared to the characteristic time evo- 
lution of the vibrons [23]. Indeed, this later time is of 
about the invert of the hopping constant $ which is gen- 
erally lesser than or about to 10 cm -1 for both molecular 



adsorbates and proteins. For instance, the dispersion of 
the stretching vibration Si-H of the system Si/H(lll) 
was found to be about 9 cm -1 [30]. In the case of the 
(2 x 1) CO monolayer adsorbed on NaCl, the lateral in- 
teraction was shown to be strongly anisotropic with a 
vibron bandwidth equal to 3 cm -1 along the longer size 
of the unit cell and to 12 cm -1 along the perpendicular 
direction [24]. In a similar way, for a- helix proteins, the 
accepted value for the vibron hopping constant is equal 
to 7.8 cm -1 (see for instance Refs. [3,19]). 

Consequently, the inequality $ x r c << 1 allows us 
two simplifying approximations in Eq.(21). First, we as- 
sume that the thermal bath does not exhibit spatial cor- 
relations ((Au n (t)Aui n >(0)) = (Au) n (t)Au) n (0))S nn i) and 
second, we neglect the non-diagonal elements of the free 
propagator exp(—iH v t). The second assumption means 
that over a time scale of about r c , a vibron does not have 
enough time to make a transition to a neighboring site. 

As a result, by inserting Eqs.(5) and (6) into Eq.(21) 
and by performing both the time integration and the 
thermal average, the TVBS relaxation rate is finally ex- 
pressed as 



dftJ(fi)[n(fi)p£ (wu(fc) + H) 



+(n(fi) + l)p£(w B (fc)-fi)] 



(22) 



where n(£Y) = (exp(—Cl/kT) — 1) _1 denotes the Bosc- 
Einstein phonon distribution at temperature T and 
where pf (w) is the projected two-vibron density of 
states defined in terms of the two-vibron Green opera- 
tor G{lo) = (ui — Hy)^ 1 as 



Pk 



(w) = -Im^- 25Z<fc',m| G(w) \ k', m) 



nN 

^ B (m)cos({k-k')m/2) 



m k 1 



(23) 



The projected density of states (DOS) p B (uj) represents 
the coupling distribution between the TVBS with wave 
vector k and energy u's(fc) and all the other two-vibron 
states which the energy ranges between ui and uj + duo. 
Such a coupling is mediated by the interaction between 
the vibrons and the phonons of the thermal bath. 

Note that the projected DOS can be expressed as the 
sum of two contributions, i.e. p B {w) = Ph^" 8 ^) + 
p B ^ F (ui), where p B ^ B (uj) and p B ^ F (uj) characterize the 
coupling distribution with the TVBS and the TVFS, re- 
spectively. The projected DOS p B ^ B (uj) can be eas- 
ily determined by considering the restriction of the two- 
vibron Green operator to the TVBS in Eq.(23). It is 
expressed as 



B^B 

Pk 



ttN ^ uj - uj B (k') + i0+ 

m k' 



$ B {m)cos{(k- k')m/2) | 2 



(24) 



Note that the knowledge of both p B (ui) and p B ^ B (uj) 
allows us to calculate p B ~* F (uj). As a result, since the 
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TVBS relaxation rate depends linearly on the projected 
DOS, it can be expressed as the sum of two contributions 
as 



■W ( k %{B^B) + W { i} 3 (B^F) (25) 



where W^(B -» B) and W^ B {B F) characterize the 
relaxation of the TVBS with wave vector k and energy 
ws(fc) over all the other bound and free states, respec- 
tively (Fig. 2b). 

In order to compare the single and two-vibron relax- 
ation rates, let us mention that the rate for a singlc- 
vibron state with wave vector q and frequency Lo(q) can 
be straightforwardly determined by using the previous 
procedure. Indeed, after performing some algebraic ma- 
nipulations, we thus obtain 

= 2n J dClJ(n)[n(Cl)g{u)(q) + fi) 

+(n(fi) + l)g(u{q) - fi)] (26) 

where g(u) stands for the single-vibron density of states 
defined as 



1 



1 



2tt$ 



"-"o ^2 
2$ > 



(27) 



As shown in Eqs.(22) and (26), the rates depend on 
the temperature via their dependence with respect to 
the average number of phonons. Moreover, both the 
two-vibron and single-vibron rates exhibit two contribu- 
tions connected to the absorption (term proportional to 
n(fi)) and to the emission (term proportional to n(fi) + l) 
of an acoustic phonon, respectively. As pointed out by 
Jakob and Persson [21], a factor 4 discriminates between 
the two rates. However, Eqs.(22) and (26) show that 
the main difference between the dephasing mechanism 
related to single and two-vibron states appears through 
the dependence of the rates with respect to the single and 
two-vibron density of states. Although the single-vibron 
density of states can be calculated exactly (Eq.(27)), the 
projected DOS p k (w) still remains unknown. Therefore, 
the next section is devoted to its calculation based on the 
determination of the two-vibron Green operator [15]. 



C. Green function calculation 

Within the intermediate basis | k,m), the Green op- 
erator G(uj) = (gj — Hy)^ 1 is block diagonal and can be 
evaluated for each wave vector k. Let us define G k as the 
two-vibron Green operator connected to the restriction 
H vk of H v to the subspace {| k, m}}, m = 0, (N— 1)/2. 
It corresponds to the Green function associated to the 
Schrodinger equation Eqs.(12)-(14). Therefore, Gk can 
be determined explicitely by taking advantage of the 
equivalence between the Schrodinger equation Eqs.(12)- 
(14) and the dynamics of a single particle [15]. Such a 



determination can be performed by applying the proce- 
dure introduced by L. Dobrzynski [31] to calculate the 
response function of super-lattices and composite mate- 
rials. 

To proceed, we first consider the tight-binding prob- 
lem on the lattice shown in Fig. lb which corresponds to 
an infinite chain. This chain is denoted the ideal lattice 
by opposition to the real lattice drawn in Fig. la. Let 
Hj^ 3 denotes the tight-binding Hamiltonian which charac- 
terizes the dynamics of this ideal lattice. The associated 
Green function, G£°(w) = (w — H£°)~ 1 , can be calculated 
analytically as 



Gf{m,m') 



(o± V^T)l m - m 'l 



(28) 



where a = (u> — 2cjo)/ ( E > fe- In Eq.(28) the sign ± is cho- 
sen in order to ensure the regular behavior of the Green 
function when | m — m' |— > oo. The Hamiltonian H v k 
can be obtained by adding an operator V k to and 
by taking the restriction to the sites which belong to the 
real lattice as 



H vk =P k (H^ + V k )P k 



(29) 



where P k = J2m=o^ 2 I ^, m){k, m | is the projector onto 
the real lattice. The goal of the operator V k , known as 
the cleavage operator [31], is first to disconnect in the 
ideal lattice the sites which do not belong to the real lat- 
tice, and then to introduce the real boundary conditions. 
The cleavage operator acts in the subspace connected to 
the interface between the real and the ideal lattice, i.e. 
m = —1,0, 1. Inside this subspace, it is straightforward 
to show that the non vanishing elements of the cleavage 
operator are expressed as 



-* fc _ 

M = i -* fc -2A (V2-l)$ fe 

o (v^-i)$ fc o 



(30) 



Finally, from Eq.(29), the Green function of the real lat- 
tice is thus expressed in terms of the Green function 
G^°((jj) of the ideal lattice, as 

G fe = P k {G? + G?V k P k (l - G^V^PkG^Pk (31) 

Therefore, by inserting Eqs. (28) and (30) into Eq.(31), 
the two-vibron Green operator can be computed exactly. 

At this step, the knowledge of the two-vibron Green 
function allows us to determine the two-vibron projected 
DOS (Eq.(23)) and then to compute the TVBS relaxation 
rate (Eq.(22)). Such a calculation is illustrated in the 
following section. 



IV. NUMERICAL RESULTS 

In this section, the previous formalism is applied 
to characterize the TVBS relaxation rate for a one- 
dimensional molecular lattice. The lateral force constant 
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$, fixed to 1.0, will be used as the frequency unit. Since 
the vibron dispersion in low-dimensional molecular lat- 
tices is rather small (lesser than or of about 10 cm -1 ), 
we assume that it can be neglected with respect to the 
phonon bandwidth, i.e. f2 c >> $. Moreover, for small 
molecules, the anharmonic parameter is usually close to 
the gas phase value and ranges between 10 — 30 cm -1 . 
We thus consider that a reasonable value for the max- 
imum of the reduced parameter A/<& is of about 10. 
Finally in order to compare the single and two-vibron 
relaxation rates, the single-vibron state with zero wave 
vector (q = 0) and with frequency u(q = 0) = ujq + 2$ 

will be used as the reference for the calculation of Wfj= 
(Eq.(26)). 

The behavior of the projected DOS p^(co) as a func- 
tion of the anharmonic parameter A is shown in Fig. 
3 for two values of the wave vector k — (full line) 
and k — ir/2 (dotted line). The projected DOS ex- 
tends from 2w - 2\M 2 + 4$ 2 to 2ui + 4$. It is max- 
imum in the frequ ency range of the bound states, i.e. 
2uj - 2y/A 2 + 4<f> 2 < u < 2ui - 2A, and exhibits two 
peaks located at the sides of the TVBS band. These 
peaks characterize the Van Hove singularities which oc- 
cur when the group velocity of the TVBS vanishes. When 
A < 2, the TVBS band and the TVFS band overlap each 
other. As a result, the projected DOS forms a single 
band which is significant over the entire frequency range. 
By contrast, when A > 2, the projected DOS exhibits 
two bands located in the frequency range of the TVBS 
and TVFS, respectively. As increasing the anharmonic- 
ity the weight of the projected DOS in the frequency 
range of the TVFS decreases drastically. Note that the 
projected DOS becomes almost independent on the value 
of the wave vector k when A » 2. 

(2) 

The behavior of the TVBS relaxation rate Wj^J B as 
a function of the temperature is shown in Fig. 4a for 
different values of the anharmonicity A = 1,3,5,7 and 
for a bath dimension a = 1. At very low temperature, 
i.e. kT « $, the rate exhibits a powerlaw dependence 
with respect to the temperature. By contrast, as increas- 
ing the temperature, the rate behaves linearly versus T 
with a slope which increases as the anharmonicity A in- 
creases. In Fig. 4b, the temperature dependence of the 
ratio Wf® B /W^ is shown for the same set of values of 
the anharmonicity. The various curves exhibit two dis- 
tinct regimes depending on the temperature. At low tem- 
perature, the ratio increases in a quasi-linear way with a 
slope which increases with A. Then, at higher tempera- 
ture, the ratio becomes almost temperature independent 
and reaches a limiting value which strongly depends on 
the anharmonic parameter. It increases as the anhar- 
monicity increases and evolves from 2.79 when A = 1.0 
to 3.85 when A = 7.0. Note that the transition between 
the two previous regimes is independent on the anhar- 
monicity and occurs when fcT/$ is of about 15-20. 

As shown in Figs. 5a and 5b, both the relaxation 

rate Wp® B and the ratio wf® B /W^ exhibit a sim- 



ilar temperature dependence when the dimension of the 
thermal bath is set to a — 2. The transition between the 
two regimes occurs at almost the same temperature, i.e. 
fcT/$ w 15 — 20 (Fig. 5b). However, as shown in Fig. 
5a, the linear evolution of the rate with respect to the 
temperature is characterized by a slope which decreases 
as the anharmonicity increases, in marked contrast with 
the previous situation. In addition, the limiting values 
reached by the ratio Wj® B /W^ at high temperature 
(Fig. 5b) decreases as A increases and evolves from 4.42 
when A = 1.0 to 1.56 when A = 7.0. Note that both 
the relaxation rate Wf® B and the ratio Wf® B /W^ 
exhibit a similar behavior when the bath dimension is set 
to a — 3 (not drawn). 

When the bath dimension is equal to a — 1, the behav- 
ior of the relaxation rate with respect to the anharmonic 
parameter A is shown in Fig. 6a for kT/& = 1 (cir- 
cle), 10 (square), 50 (triangle), 100 (diamond) and for 
k = (full line) and k — n/2 (dotted line). Whatever 
both the temperature and the value of the wave vector, 
the ratio W^ B /W { q _q increases with the anharmonicity 
and reaches a limiting value for a strong anharmonicity. 
At low temperature, i.e. kT = $, this limiting value 
is almost equal to unity whereas it appears close to 4 
at high temperature. Note that the rate for the TVBS 
with wave vector k = is lesser than the rate for the 
TVBS with wave vector k = n/2. Nevertheless, the dif- 
ference is less and less important as the anharmonicity 
becomes stronger. In Fig. 6b, the behavior of the ratio 

W k%( B -> F )/ W k%( B B ) with rcs P ect to A is shown 
for the same set of temperatures. For both k = and 
k = n/2, the ratio decreases as the anharmonicity in- 
creases and increases with the temperature. It is almost 
lesser than unity excepted at high temperature and low 
anharmonicity. Note that when kT /<& > 10, the curves 
decrease with respect to A according to an invert power- 
law. For instance, when kT = 50$ and k — 0, the ratio 
evolves from 1.02 when A = 1 to 0.014 when A = 10. 

As shown in Figs. 7a and 7b, the relaxation rate be- 
haves in a fully different way when the bath dimension is 
set to a — 2. Indeed, although the ratio wf^g/wfj= ap- 
pears almost independent on the anharmonicity at very 
low temperature (kT = $), it decreases as the anhar- 
monic parameter A increases when fcT/$ = 10, 50, 100 
(Fig. 7a), in a marked contrast with the previous results. 
For instance, when kT = 50$ and k = 0, the TVBS re- 
laxation rate is greater than the single-vibron relaxation 
rate by a factor 4.17 when A = 1 whereas this ratio de- 
creases to 0.98 when A = 10. In addition, the relaxation 
rate for the TVBS with wave vector k = appears now 
greater than the relaxation rate for the TVBS with wave 
vector k = tt/2 but, as previously, the difference is less 
and less important as the anharmonicity increases. As 
shown in Fig. 7b, the ratio wf ] B {B -» F)/W^ B {B -» B) 
decreases as the anharmonicity increases. However, it is 
greater than unity over a large anharmonicity range, ex- 
cepted at low temperature. For instance, when kT = 10$ 
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and k = 0, the ratio, which is equal to 3.31 when A = 1, 
becomes lesser than unity when A > 4.0. Note that, 
the various curves exhibit two regimes depending on the 
anharmonicity since the ratio rapidly decreases at low 
anharmonicity, i.e. A < 2.5, whereas it decreases more 
slowly at high anharmonicity. 

The same results are shown in Figs. 8a and 8b when 
the bath dimension is equal to a — 3. As in the previous 
case, the ratio W^/W^ is almost independent on the 
anharmonicity at very low temperature (kT = $), but 
decreases as the anharmonicity increases when fcT/$ = 
10,50,100 (Fig. 8a). For instance, when kT = 50$ 
and k = 0, the TVBS relaxation rate is equal to 7.01 
time the single-vibron rate when A = 1 whereas this 
ratio is equal to 3.89 when A = 10. Moreover, the re- 
laxation rate for the TVBS with wave vector k = is 
greater than the relaxation rate for the TVBS with wave 
vector k = ir/2. Note that, at high temperature, the 
rate clearly shows two distinct regimes since it exhibits 
a fast decrease at low anharmonicity (A < 3.5) whereas 
it decreases more slowly at high anharmonicity and be- 
comes almost constant. As shown in Fig. 8b, the ratio 

W i%( B -* F )/ w t%( B -» B ) behaves in a fully differ- 
ent way with respect to the two previous situations. At 
low temperature, i.e. kT = $, it is lesser than unity and 
decreases as the anharmonicity increases. When the tem- 
perature is equal to kT = 10$, it first slightly decreases 
at very low anharmonicity (A < 2), then it increases with 
A to reach a maximum value equal to 10.45 at A w 7.0. 
Finally, as increasing the anharmonicity from the previ- 
ous value, the ratio decreases. At higher temperature, 
i.e. fcT/$ = 50 and 100, the ratio shows a similar be- 
havior but increases more drastically. Note that it does 
not reach a maximum value into the anharmonicity range 
considered in this work. 



V. DISCUSSION 

To interpret and discuss the previous numerical results, 
let us first focus our attention on the temperature de- 
pendence of the TVBS relaxation rate. Indeed, what- 
ever both the anharmonicity and the bath dimension, 
the rate exhibits a powerlaw and a linear dependence on 
the temperature at low and high temperature, respec- 
tively. These features are due to the fact that the relax- 
ation involves single phonon processes in the course of 
which a phonon is absorbed or emitted. As a result, the 
rate depends linearly on the phonon Bosc-Einstcin dis- 
tribution which is responsible for the observed behavior. 
When k — 0, the TVBS frequency lies at the bottom of 
the TVBS band. The relaxation involves phonon absorp- 
tion only, and, therefore, the corresponding rate vanishes 
as the temperature tends to zero since the number of 
phonons becomes infinitcsimally small. By contrast, for 
a finite wave vector, phonon emission can occur leading 
to a finite value of the rate at zero temperature. Note 



that in our simulation we used a positive value for the 
lateral force constant $. As a result, the single vibron 
state q = 0, which was considered as the reference for the 
calculation of the rate W^Iq, lies at the top of the single- 
vibron band. Therefore, W^} reaches a finite value at 
T = since phonon emission can takes place and depends 
linearly on the temperature at high temperature. 

The more surprising results are clearly related to the 
behavior of the TVBS relaxation rate with respect to 
the anharmonicity. Indeed, whatever the value of the 
bath dimension, the rate is a slowly varying function of 
the anharmonicity at very low temperature, and the re- 
laxation path involves the decay of the TVBS over the 
other TVBS, only (See Figs. 7b, 8b and 9b, for kT = $). 
However, at sufficiently high temperature (kT > 10$), 
the numerical results indicate that the rate exhibits two 
distinct regimes depending on the bath dimension. When 
a = 1, the TVBS relaxation rate increases with the an- 
harmonicity and reaches a limiting value almost four time 
greater than the single-vibron rate. Moreover, the main 
contribution of the rate involves the relaxation over the 
other TVBS, only. In a marked contrast, when the bath 
dimension is equal to 2 and 3, the rate decreases as the 
anharmonic parameter A increases. When a = 2, the re- 
laxation involves the decay of the TVBS into both TVBS 
and TVFS, the two processes contributing in a similar 
way to the rate. By contrast, when a = 3, the decay 
into the TVBS becomes negligible and the main contri- 
bution of the rate involves transitions from the TVBS to 
the continuum of TVFS. 

Such features originate in the convolution product be- 
tween the projected DOS p^f(w) and the coupling dis- 
tribution J(ui) weighted by the Bose-Einstein distribu- 
tion. When kT « $, the Bose-Einstein distribution 
decreases exponentially with the phonon frequency and 
thus favors the exchange of very low frequency phonons. 
Therefore, whatever the bath dimension, the TVBS can- 
not decay into the TVFS continuum and relaxes over the 
other TVBS. When kT >> $, both phonon emission and 
phonon absorption contribute similarly to the rate since 
n(f2) i=a n(tt) + 1 kT/VL. The rate reduces to the convo- 
lution between the projected DOS and the distribution 
kTJ(fl)/n. This latter function scales as fi Q_1 and is 
constant when a = 1 whereas it increases with fi when 
both a = 2 and 3. 

Therefore, when a = 1, the relaxation of the TVBS 
mediated by phonon exchange occurs with the same 
probability whatever the value of the frequency of the 
phonon involved in the process. As a consequence, the 
behavior of the rate is governed by the strength of the 
coupling between the TVBS and the other two-vibron 
states induced by the vibron/phonon interaction. From 
Eqs.(22)-(23), this coupling, which is expressed as the 
integration of the projected DOS, leads to a relaxation 
rate written as 
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<1 



2irkTJ 



Y,\^k(m)co S ((k-k')m/2)\ 2 (32) 



After some straightforward calculations, Eq.(32) leads to 



a ratio wj. 2 B /Wg= expressed as 



(i) 



W, 



(2) 



k,B 



2A 



W, 



(i) 



9=0 



y/A 2 + (2$cos(fc/2)) 2 



(33) 



The ratio W^ 2 B /W^ increases with the anharmonicity, 
in perfect agreement with the numerical results shown in 
Fig. 6a , and evolves from 2 when A = to 4 when A 
tends to infinity. In other words, for a small anharmonic- 
ity, the TVBS relaxation rate is twice the single- vibron 
rate which indicates that the two vibrons are almost in- 
dependent, i.e. the TVBS rate appears as the sum of the 
rates connected to each vibron. As increasing the an- 
harmonic parameter, the length of the bond between the 
two vibrons (Eq. (17)) reduces and the local nature of 
the vibron population increases to reach a limiting value 
equal to 2. The strength of the coupling, proportional to 
the square of the population, leads to a factor 4 between 
the TVBS rate and the single-vibron rate. Note that as 
increasing the value of the wave vector, the local nature 
of the TVBS increases since the length of the bond be- 
tween the two vibrons decreases (cf Eq.(17)) leading to 
an enhancement of the relaxation rate. 

In a marked contrast, when a = 2 and 3, the behavior 
of the TVBS relaxation rate results in the competition 
between the projected DOS and the coupling distribu- 
tion kTJ(fl)/fl which contribute in an opposite way. In- 
deed, the projected DOS, which is maximum in the fre- 
quency range of the bound states, favors transitions from 
the TVBS to the other TVBS. Such transitions can thus 
be induced by the exchange of low frequency phonons. 
However, the coupling distribution kTJ(il)/il, which in- 
creases with the phonon frequency, favors transitions me- 
diated by high frequency phonons from the TVBS to the 
TVFS continuum. At high temperature, the numerical 
results clearly show that the decay of the TVBS into the 
TVFS continuum is the main contribution of the rate. 
This feature indicates that the influence of the phonon 
coupling distribution prevails. However, the contribution 
of the projected DOS in the frequency range of the TVFS 
decreases drastically as the anharmonicity increases and 
counterbalances the growth of the phonon coupling dis- 
tribution. As a result, the rate decreases with respect 
to the anharmonic parameter. The effect of the phonon 
coupling distribution becomes stronger as the bath di- 
mension increases. Indeed, when a — 2, both decays into 
the TVBS and TVFS continuum, occur with almost the 
same probability. However, when a = 3, the latter path 
for the relaxation becomes strongly dominant. Note that 
as increasing the value of the wave vector, the energy 
difference between the TVBS and the TVFS continuum 
decreases. As a result, the frequency of the phonon ex- 
changed during the relaxation mechanism decreases lead- 



ing to a lesser value for the relaxation rate since the cou- 
pling distribution J(f2) is reduced. 

When the anharmonicity is sufficiently strong, an an- 
alytic expression of the projected DOS in the frequency 
range of the TVBS can be determined. Indeed, by per- 
forming a second order expansion of Eq.(16) with re- 
spect to Q/A, the frequency of the TVBS is expressed 
as u B (k) w 2uj -2A - 2($ 2 /A)(l + cos(k)). There- 
fore, the TVBS dynamics can be compared with the dy- 
namics of a single vibron characterized by a force con- 
stant = -$ 2 /A The projected DOS p B ^ B (uj) is thus 
given by the single vibron density of states Eq.(27) ac- 
cording to the correspondences <& — > <P 2 /A and cj — > 
2lu -2A- 2<5> 2 /A. The rate W^ ] B (B B) for the de- 
cay of the TVBS into the other TVBS is thus given by 
Eq.(26) with the previous correspondences and by adding 
a factor four to account for the population effect. For in- 
stance, at high temperature, the rate for the relaxation 
of the TVBS with zero wave vector is written as 



where 1(a) is defined in terms of the Gamma func- 
tion as 1(a) = 2 3a T(a + l/2) 2 /T(2a + 1). As shown 

(2) 

in Eq.(34), the rate Wj^J B (B — > B), which scales as 
A 1 ^ 01 , is constant when a = 1 whereas it decreases when 
a > 1. Although no analytical expression for the rate 
Wf, B (B — > F) can be determined, the numerical results 
clearly show that this rate scales as 1/A when a = 1 and 
becomes negligible compared with wjf B (B — > B). By 

contrast, when a = 2 and 3, the rate W^ B (B — > F) de- 
creases with A according to a law decaying more slowly 
than A 1 ^", and is the dominant contribution of the rate. 

At this step, let us mention that the numerical results 
presented in the previous section were obtained by con- 
sidering integer values for the bath dimension. However, 
the formalism can be applied to real values of the bath di- 
mension. In particular, at high temperature and strong 
anharmonicity, our calculations predict that the TVBS 
relaxation rate decreases with the anharmonicity when 
a > 1 . In addition, the dominant relaxation pathway cor- 
responds to the decay of the TVBS into the other TVBS 
over a large a range, i.e. 1 < a < a c where a c < 2. In 
that case, the total relaxation rate is thus approximately 
given by Eq.(34). For instance, when kT = 100$ and 
A = 10, the ratio W^ ] B (B -» B)/W^ 2) B (B -» F) is equal 
to unity for a bath dimension a c w 1.9. 

To illustrate the previous theoretical calculations, let 
us compare our results with experimental data. The sys- 
tem CO/NaCl(100) was previously studied using molec- 
ular dynamics simulations [24]. The CO monolayer has 
a (2 x 1) structure with a rectangular unit cell contain- 
ing two molecules mutually perpendicular. Since the lat- 
eral interaction was shown to be strongly anisotropic, 
the vibron dynamics of the monolayer exhibits a one- 
dimensional character. For instance, along the longer size 
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of the unit cell, the lateral force constant between near- 
est neighbor molecules is positive and equal to <f> = 2.94 
cm -1 . Assuming an anharmonic parameter equal to the 
gas phase value, i.e. A = 13.3 cm -1 , the reduced an- 
harmonicity is A/<I> = 4.5. At T = 5 K and for a bath 

dimension a = 1, the theoretical ratio wj^} Q B /W^ is 
equal to 1.3. Note that for a bath dimension equal to 2 or 
3, the ratio is drastically reduced by a factor of about 10. 
This result is in a perfect agreement with the experimen- 
tal ratio, equal to 1.55, obtained by Dai and Ewing [11]. 
A genuine one-dimensional lattice was obtained by Jakob 
[14] by carefully controlling the adsorption of CO on a 
Pt(lll) surface exhibiting steps. As measured by the 
author, the dynamical coupling between the molecules is 
responsible for a blueshift Alu = 19 cm -1 . By assum- 
ing nearest neighbor interactions, this coupling leads to 
a positive lateral force constant <f> = 9.5 cm -1 . In ad- 
dition, the anharmonicity of the system, measured by 
Jakob, yields a reduced parameter A/& = 1.42. By fol- 
lowing the experiment performed by the author, we thus 
obtain a ratio Wj® B /W^ = 2 at T = 80 K, which is 
in good agreement with the experimental ratio equal to 
1.5. 

To conclude, let us discuss the main implications of 
the previous results on the energy transport mediated by 
TVBS. Indeed, from a quantum kinetic point of view, the 
vibron/phonon interaction is responsible for the breaking 
of the coherence and favors a diffusional motion of the 
vibrons after a time scale of about the dephasing time 
[23]. The previous results clearly show that such a dif- 
fusional motion, for TVBS, drastically depends on the 
bath dimension. When a = 1, the coupling with the 
phonons is not strong enough to break the bond between 
the two vibrons which cannot decay into the continuum 
of TVFS. Therefore, the thermal bath keeps the internal 
coherence between the two vibrons but induces momen- 
tum scattering which is responsible for transitions from a 
TVBS to another TVBS. Such a scattering in momentum 
space corresponds to incoherent hops in real space in the 
course of which the center of mass of the two vibrons has 
a diffusional motion. Therefore, although the energy is 
transported in an incoherent way, it diffuses by packets 
along the lattice. Note that since the TVBS relaxation 
rate is almost four time greater than the corresponding 
single- vibron rate, a similar ratio is expected for the dif- 
fusion coefficients. By contrast, when a — 2 or 3, the 
phonon bath favors the breaking of the bond between 
the two vibrons over a large anharmonicity range. As 
a result, the decay of the TVBS into the continuum of 
TVFS takes place and the two particles will tend to dif- 
fuse independently leading to the uniformization of the 
energy along the lattice. 

Finally, let us mention some remarks on forthcoming 
work. First, the present results were obtained by con- 
sidering the coupling between the vibrons and acoustic 
phonons. However, as discussed in sec. II, the formal- 
ism can be applied to different kinds of thermal bath. 



In particular, it would be interesting to investigate the 
effect of optical phonons which can be responsible for 
the occurrence of resonances between the TVBS band 
and the TVFS continuum. These features were pointed 
out experimentally by Jakob et al. [21] when studying 
the linewidth of the overtone of the system CO /Ru and 
attributed to the coupling with the frustrated transla- 
tion mode of the monolayer. Second, we have restricted 
our analysis to processes involving a single phonon, only. 
However, two-phonon processes can contribute signifi- 
cantly to the relaxation rate, depending on the bath di- 
mension. For instance, two-phonon processes may en- 
hance the spectral distribution of the vibron/phonon cou- 
pling at low frequency favoring the decay of TVBS into 
the TVBS band. 
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Figure Caption 

Figure 1 : Equivalence between the Schrodinger equa- 
tion of two-vibron states in the intermediate basis and 
the dynamics of a single fictitious particle moving quan- 
tum mechanically in a one-dimensional lattice (see the 
text). The index m, which characterize the separating 
distance between the two vibrons, labels the sites of the 
equivalent lattice, (a) The equivalent lattice, called the 
real lattice, extends from m — to m = (N — l)/2. (b) 
Ideal lattice used to calculate the two-vibron Green op- 
erator. The circles define the self-frequencies of the sites 
m and the lines characterize the hopping constants 
(thin line) and V%&k (full line) (see the text). 

Figure 2 : (a) Two-vibron eigenstates energy spectrum 
which exhibits two bands. The high frequency band, 
located around 2w 0j contains the two-vibron free states 
continuum. The low frequency band is the dispersion re- 
lation of the two-vibron bound states. The anharmonic 
parameter is equal to A = 2.5$. (b) Schematic represen- 
tation of the relaxation processes. 

Figure 3 : Projected DOS (oj) as a function of the 
anharmonic parameter A for two values of the wave vec- 
tor k = (full line) and k — ir/2 (dotted line). 

Figure 4 : (a) Temperature dependence of the two- 

vibron bound state relaxation rate W^_ B for different 
values of the anharmonicity A = 1 (circle), 3 (square), 5 
(trinagle), 7 (diamond) and for a bath dimension a = 1. 
(b) Temperature dependence of the ratio Wj^ B /W^ 
for different values of the anharmonicity A = 1 (circle) , 
3 (square), 5 (trinagle), 7 (diamond) and for a bath di- 
mension a = 1. 

Figure 5 : (a) Temperature dependence of the two- 

vibron bound state relaxation rate W^J B for different 
values of the anharmonicity A = 1 (circle), 3 (square), 5 
(trinagle), 7 (diamond) and for a bath dimension a = 2. 
(b) Temperature dependence of the ratio wj^} B /W^ 
for different values of the anharmonicity A = 1 (circle), 
3 (square), 5 (trinagle), 7 (diamond) and for a bath di- 
mension a = 2. . 

Figure 6 : (a) Behavior of the TVBS relaxation rate 
with respect to the anharmonicity A for fcT/$ = 1 (cir- 
cle), 10 (square), 50 (triangle), 100 (diamond) and for 
k = (full line) and k = tt/2 (dotted line). The bath 
dimension is equal to a = 1. (b) Behavior of the ratio 

W k%( B -> F )/ W k%( B -> B ) with respect to the an- 
harmonicity A for the same set of parameters as in Fig 
6a. 

Figure 7 : (a) Behavior of the TVBS relaxation rate 
with respect to the anharmonicity A for kT / $ = 1 (cir- 
cle), 10 (square), 50 (triangle), 100 (diamond) and for 
k = (full line) and k = tt/2 (dotted line). The bath 
dimension is equal to a = 2. (b) Behavior of the ratio 

W k%( B -> F )/ W k%( B -> B ) with aspect to the an- 
harmonicity A for the same set of parameters as in Fig 
7a. 



Figure 8 : (a) Behavior of the TVBS relaxation rate 
with respect to the anharmonicity A for kT j $ = 1 (cir- 
cle), 10 (square), 50 (triangle), 100 (diamond) and for 
k = (full line) and k = tt/2 (dotted line). The bath 
dimension is equal to a = 3. (b) Behavior of the ratio 

W kl( B -» F )/wi%(B -» B) with respect to the an- 
harmonicity A for the same set of parameters as in Fig 
8a. 
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